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Abstract. Here we report on the first successful exoplanet transit observation with the Stratospheric Observatory for 
Infrared Astronomy (SOFIA). We observed a single transit of the hot Jupiter HD 189733 b, obtaining two simultaneous 
primary transit lightcurves in the B and z' bands as a demonstration of SOFIA’S capability to perform absolute transit 
photometry. We present a detailed description of our data reduction, in particular the correlation of photometric 
systematics with various in-flight parameters unique to the airborne observing environment. The derived transit depths 
at B and z' wavelengths confirm a previously reported slope in the optical transmission spectrum of HD 189733 b. Our 
results give new insights to the current discussion about the source of this Rayleigh scattering in the upper atmosphere 
and the question of fixed limb darkening coefficients in fitting routines. 
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1 Introduction 


The exploration and detailed analysis of exoplanet atmospheres is one of the most dynamic fields of 
astrophysics today. The first successful observations probing the atmospheri^roperties of Jupiters 
were conducted with the Hubble Space Telescope’s (HST) STIS instrument (UJ) and with the Spitzer 
Space Telescope’s (Spitzer) IRS instrument (Pilsl). Approximately a decade later, we are now able 
to analyze the atmospheres of planets down to Neptune (GJ 436b, Idj) and super-Earth (GJ 1214b, 
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51) sizes. High signal-to-noise observations of the absorption depth of moleeular bands with strong 
signatures in the optieal and infrared portions of the speetrum ean determine both the temperature- 
pressure thermal profile as well as the abundanees of atmospherie eonstituents, helping to eonstrain 
the overall atmospherie ehemistry and strueture ([y; lZ|; IS]). It has also been shown that the relative 


abundanees of atomie speeies sueh as C, O, and N in a gas giant’s atmosp 


lere eould be indieative 


of the region of the proto-planetary disk in whieh the planet formed Ml id: I 111). 

In this paper we present the first observations of a transiting planet, the well-known transiting 
Hot Jupiter HD 189733 b, using the Stratospherie Observatory for Infrared Astronomy (SOFIA, 
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131). SOFIA eonsists of a 2.5-meter teleseope mounted within a modified Boeing 747-SP air- 


eraft operating at altitudes up to 45,000 feet, thereby offering the opportunity for observations at 
altitudes where the tellurie absorption from the Earth’s atmosphere is greatly redueed. In Seetion[2] 
we deseribe our observations, in Seetion[3]we deseribe our data reduetion and transit light eurve fit¬ 
ting analysis, in Seotion|4]we deseribe our results, and diseuss the implieations for the atmosphere 
of HD 189733 b and for future observations with SOFIA in Seotion[5l 


1.1 Exoplanet observations with SOFIA 

As an airborne observatory, SOFIA has a number of potential advantages for preeise time-domain 
speetrophotometrie observations of transiting exoplanets. Ground-based observations are signif- 
ieantly affeeted by variations of absorption from tellurie gases, in partieular H 2 O, in the Earth’s 
atmosphere, and these same gases are also the speeies of interest in exoplanet atmospheres' SOEIA 
ean observe in important atmospherie windows not observable from the ground (Il4l: Il5l: Il6l: IITI) . 
These are mostly the water bands but also CO, CH 4 , CO 2 are mueh better mixed and therefore 


reduee the temporal variation in these bands, whieh is a erueial point for time-series observations. 
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These bands are also the molecular bands of interest in targets such as HD 189733 b. The air¬ 


borne observatory operates in the wavelength regime where the planet’s black-body temperature 
peaks and contrast ratios between star and planet improve, and the SOFIA telescope also operates 
at lower temperatures (240K) than ground-based telescopes; therefore the contribution from the 
thermal background (the dominant noise source for transit observations at wavelengths longer than 
3 microns) are significantly reduced. SOFIA can observe simultaneously at infrared and optical 


wavelengths using its FLITECAM (1181) and HIPO (1191) instruments in ’FLIPO’ mode (l20h . thereby 


obtaining light curves for a single transit event over a wide range in wavelength. 

However, there are also certain challenges when observing with SOFIA due to the airborne 
observing environment. In test flights HIPO photometry has been shown to be affected by terres¬ 
trial Rayleigh scattering, by ozone extinction in the Chappuis band that varies with position along 


the flight 


path, potentially by volcanic aerosols, by pointing errors, and by other factors affecting 


the PSF (1211) . At optical wavelengths the PSF is dominated by wavefront aberrations imposed 


by the turbulent shear layer that passes over the telescope cavity. The strength of the shear layer 
density fluctuations varies with the static air density, causing density-dependent interaction of the 
very broad wings of the PSF with the photometric aperture. Furthermore, the logistics of aircraft 
operations result in observation windows which are limited in duration and dependent on the flight 
plan chosen for a specific flight; this can result in the limited availability of measurements before 
or after a transit event and hamper the decorrelation of these various observational effects with the 
transit light curve. 

Fortunately an expected dependence on Mach number has not been seen, and the residual 
systematic noise due to impacts of focus errors and high speed image jitter are much less than the 
photon noise if large circular synthetic photometric apertures are used (see also discussion in 14.II) . 
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Other potential sourees of systematic error that have not yet been investigated include higher-order 
extinction corrections, water vapor absorption in the z' band, possible polarization effects from 
the SOFIA tertiary mirror, and the impact of a known temperature-dependent astigmatism term 
in the PSF. We have also not yet investigated whether weighted aperture photometry is beneficial 
for precise photometric work. Furthermore, Rayleigh extinction is well behaved in flight because 
the aircraft flies at constant pressure (and therefore constant Rayleigh zenith optical depth) and the 
usual extinction correction can be applied. Ozone extinction can largely be avoided by proper filter 
selection and volcanic aerosols are thankfully uncommon. Pointing and focus errors, and variations 
with instrument and environmental factors must be decorrelated using SOFIA housekeeping data, 
to reach the high signal to noise needed for our differential transit depth measurements. 


HD 189733 b 


HD 189733 b (l22h is a Jupiter-mass planet orbiting at 0.03 AU around one of the closest K-type 


stars; the deep transit signal in addition to a very bright, nearby host star results in the best op¬ 
portunity for high-precision characterization of any known exoplanet. This transiting system is a 
benchmark for exoplanet observations and has been the target for many ground- and space-based 
observations. Both multi-band photometry and spectroscopy with Spitzer have provided measure¬ 
ments of the mid-IR emission from the planet by measuring the occultation of the planet by the 
central star, probing both molecular absorption and the temperature structure in the bulk of the up¬ 
per atmosphere (1231: l24t) . Observations with HST have explored molecular bands in the NIR, and 
indeed, early results with the NICMOS spectrograph claimed absorption from H 2 O (l25l:l26h. (3*774 
(l27h and CO 2 between 1.5 and 2.5 /im . Other observations found the same molecular features 
in the emission spectra observed during secondary eclipse observations dj; 1281) . However, there 
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has been eonsiderable diseussion on this topie and the quality and reproduetion of these results 


(e.g., l29|; l30|; l3lh . It is still uneertain whether these moleeular features exist or if there is a haze 


that obseures wavelengths below 2 /rm (1321:13311. in faet, (13411 make the argument that HD 189733 


b’s atmos 


rhere is most likely dominated by Rayleigh seattering in the visible and near-infrared. 


While (l35b an d other observations did not direetly deteet any starspots by erossing events, a reeent 


study by (12611 argued that the measured slope at shorter wavelength eould also partly be eaused by 


unoeeulted star spots, redueing the eontribution of Rayleigh seattering to moleeular hydrogen and 
not neeessarily dust in HD 189733 b’s atmosphere. 

Our SOFIA observations were designed to both examine the presenee of a strong Rayleigh 
slope in the optieal (using photometrie observations in B and z' filters with HIPO) and either eon- 
firm or rejeet the existenee of absorption from H 2 O in the NIR (using photometrie measurements 
in FLITECAM’s Pasehen alpha 1.88 /im filter, just longwards of the eurrent upper limit of HST 
and unobservable from the ground). While both of these features were already measured in various 
observations with HST WFC3 and NICMOS (see referenees above) they still added a eompelling 
seientifie value (reproduetion and direet eomparison) to our observations, whieh were mainly a 
proof of eoneept experiment for SOFIA. These observations leverage the advantages of SOFIA 
for simultaneous optieal and near-infrared observations that are diffieult or impossible from the 
ground, as well as providing an optimal target for initial tests of the preeision possible for exo¬ 
planet transits with SOFIA. 


2 Observation 

We observed HD 189733 b during a transit on SOFIA’S flight number 134 on UT Oet 1 2013 as 
part of a Cyele 1 GO program (PI: Mandell, Proposal ID: 01-0099). Observations were eondueted 


5 



















in the FLIPO configuration (FLITECAM and HIPO operating simultaneously) in order to observe 
in three optical and infrared bands at the same time: B and z' with HIPO and a narrow-band 
filter covering the Paschen a spectral feature at 1.88 /rm with FLITECAM. The HIPO filters were 
selected to avoid spectral regions with potentially high ozone variability, while the FLITECAM 
filter was chosen due to its wavelength coverage of a prominent H 2 O spectral feature that cannot be 
sampled from ground-based observatories. An additional optical channel (for general calibration or 
tracing of a specific telluric absorption band) can be obtained from the Local Plane Imager (PPI- 1 -), 
though we were unable to acquire EPI-i- data for our flight due to an instrument malfunction (see 
below). 

Due to constraints on the flight plan imposed by requirements on the direction and timing of 
SOEIA flights, the observing period only allowed for a very short baseline before ingress and 
almost no baseline after egress. Two different instrument operators were running the blue and red 
sides and due to some problems with the general acquisition of stable condition for observing (see 
below), the blue side started later than the red side. In detail we got 0.00 hours before, 0.20 hours 
after, and 1.77 hours of in-transit for the blue side and 0.26 hours before, 0.20 hours after, and 
1.76 hours of in-transit for the red side. The duty cycle for both HIPO channels is essentially 100 
percent with a minimal dead time due to the CCD frame transfer of a few milliseconds, because in 
the frame transfer CCD the readout is overlapped with the next integration. 

Various instrumental issues emerged before and during the flight that further hampered the ac¬ 
quisition of high-quality data. In mid September 2013, several components in the FLITECAM 
detector computer failed and were replaced by spares. At the same time, a required change to the 
liquid helium venting configuration led to the creation of a thermo-acoustic oscillation that reduced 
the liquid He hold time by approximately 50%. Llights were modified to account for the shorter 
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Helium hold time, but unfortunately, the replaeement eomputer eomponents led to improper ini¬ 
tialization of the FLITECAM deteetor. In the ease of observations in whieh we simply stare at the 
souree, or observation sequenees utilizing eoadds, this quiekly resulted in souree saturation. The 
problem was not notieed until the eonelusion of the September/Oetober 2013 flight series and the 
start of data analysis in mid Oetober 2013. The Helium hold time as well as the eleetronies prob¬ 
lem have been fixed in the meantime, and we sueeessfully demonstrated FLITECAM observations 


of bright targets during FLIPO eommissioning and early seienee in Feb 2014 (l36t) . 


In order to eapture the light that is passed through the SOFIA teleseope’s diehroie tertiary 
mirror (25% and 45% refleetivity for the B and z’ bandpasses; seel20h to the Foeal Plane Imager 
(FPI-I-) we had planned to save its image data for seientifie analysis. Most of the visual light passes 
the tertiary beam splitter before it is refleeted into the Nasmyth tube by the fully-refleetive tertiary. 
The peak transmission of the tertiary beam splitter is at 570 nm with 64% of the light passing the 
mirror. A signifieant amount of visual light is not transmitted but rather absorbed or refleeted along 
with the longer, infrared wavelengths. However, in the range between 480 nm to 800 nm where the 
visual-light CCD eameras are most sensitive, more than 50% of the light is transmitted. The FPI-i- 
eontains a highly sensitive and fast EM CCD eamera. Its images are primarily used for traeking but 
ean also be stored without disrupting the traeking proeess and in parallel with measurements of the 
instruments mounted to the teleseope. Unfortunately, the EPI-i- suffered a failure of its eontroller 
eleetronies right at the beginning of the exoplanet observation leg, foreing us to abandon EPI-i- data 
taking and traeking. 

A faulty power supply in the eontroller eleetronies was found to be the eause of the failure. 
Sinee its replaeement the EPI-i- has worked reliably during future flights, ineluding a mission for 
an exoplanet transit observation in early 2014 (PL Angerhausen & Dreyer, in prep.), where we 
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were able to eharaeterize its photometrie performanee. Starting with SOFIA observing eyele 4 
(2015), the FPI+ will be made available for proposals as a faeility seienee instrument. 

The HIPO instrument was operated in Basie Oeeultation mode with full-frame read-out to max¬ 
imize our field standard possibilities, though we later eoneluded that the available field standards 
were too faint to use for differential photometry measurements. Our final HIPO data set eonsists 
of about ~2700 measurements of 3 s integrations on red side in the Sloan z' filter and ~1000 
measurements of 7 s integrations on the blue side in the Johnson B filter. As a baek-up for the 
FPI, SOFIA’S Fine Field Imager (FFI) was used for traeking. Due to its lower spatial resolution 
and defoeussing during the eool-down phase, traeking aeeuraey was redueed. The time required to 
reeover from the FPI-i- failure and switch to the FFI also caused nearly all of the pre-transit baseline 
to be lost. This loss was particularly acute because the flight plan, scheduled late in the observing 
season, included no post-transit baseline, leaving us in the difficult position of analyzing a transit 
with no out-of-transit baseline for noise characterization. 

3 Data reduction 

3.1 Image processing 

Due to the failure of FLITECAM and FPI only the HIPO data was reduced. The raw HIPO CCD 
frames were bias-subtracted and overscan-corrected. Flat fielding was complicated by the fact that 
the shape of twilight sky flats did not match the behavior of moving star images. We suspect that 
this was due to a stray light issue, possibly due to insufficient baffling or due to the reflective con¬ 
ical “button” installed in the central obscuration of the SOFIA telescope to reduce FLITECAM’s 
thermal background. In the B filter the 1-2% amplitude laser annealing pattern common in e2v 
CCDs was apparent, and in both filters a pattern similar to vignetting was seen in the sky flats 
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that was not seen in the positional dependenee of star flux. Our final flat fielding approaeh for 
the B filter data was to fit the overall eurvature with a low-order polynomial and divide it in the 
twilight flats leaving the laser annealing pattern intaet but removing the overall eurvature.Flat field 
eorreetion was done with this modified flat. The removal of the low order eurvature was neeessary 
to isolate the low-level annealing patterns. Otherwise the flatfielding will be dominated by the 
high-flux sky flat and introduee a bias. The flatfielding then removes the annealing pattern. On 
the red side we bypassed flat field eorreetion altogether. We are eurrently working on lab tests to 
understand and improve the flat fielding proeedure for future HIPO observations. 

3.2 Aperture photometry 

We performed aperture photometry using large aperture diameters (60 aresee) to reduce the effects 
of focus, jitter, and tracking errors. For our very bright target, sky and read noise contributions were 
each only ~ 10% of the stellar photon noise even in this very large aperture (see also section |4T]) . 
We avoided PSF fitting because the shear layer PSF is not well represented analytically. However 
we did use PSF fitting to estimate the brightness ratio of the target and its faint late-type companion 
and to derive the characteristics of the PSF shape (i.e. FWHM, ellipticity) for each exposure. The 
total number of electrons per exposure is 5.81 x 10® in the B channel, and 1.44 x 10^ in z'. The 
total number of sky electrons in this aperture is 1.46 x 10® forB, and 2.83 X 10® in x'. See section 
l4.1l for a detailed noise budget. 

The initially derived transit depths were diluted due to the inclusion of the fainter companion 
star (HD 189733 B, M4V, d~ 12”, V~ lOma^f) in our photometric aperture. The dilution was 
undetectable in the B filter and ~ 1.5% in the z' filter. Propagating this error into the final radius 
ratio, the correction for Rp/R^ should be 0.0012 at z' (see Tabled footnotes), any other stars in 
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the aperture are unmeasurable in z’. In B neither eompanion is measurable. 


3.3 Decorrelation and transit fitting 

The absolute raw lighteurves obtained with aperture photometry elearly show a transit signature 
(see Figure [I])- Correction of flux variations using the nearby field stars was not necessary. How¬ 
ever, the time series of the comparison stars were added to the base parameter set for the following 
analysis in order to allow for this correction if significant on a lower level. Table [U lists the cor¬ 
relation of all parameters and shows that the effect of the comparison stars is negligible. Other 
observational parameters such as the focus, overscan or air density have a much larger correlation: 
An initial inspection of the light curves revealed a significant amount of this correlated (’red’) noise 
in the raw photometry, such as a short increase in flux during the transit, with a shape similar to 


the deviation that would be caused by the planet transiting a cooler starspot (137 


381). To avoid 


bias from this feature in our initial light curve fitting, we fit the light curves with the signature 
masked out. Absent out-of-transit baseline data we attempted to correct for the correlated noise 


by combining an analytic model for the exoplanet transit light curve (l39t) with a linear combina¬ 
tion of normalized (through division by their mean) observational parameters pfit), sampled at the 
same time as our exposures, as a model of the residuals (Rmodeiit) = X] c* x Pi{t))- For the light 


curve model, we fixed the planetary orbital parameters to values from previous observations (140 . 


see Tables I2l and 13]) but left the transit depth and stellar limb darkening (optional) as free param¬ 
eters. The model parameters were optimized through a multiple linear regression analysis, and 
this noise model was then subtracted from the original data. Table [H shows the (linear Pearson) 
correlation coefficients of each observational parameter with the respective residuals as well as the 


mean values that were used for normalization. 
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HlPO B extinc. ond flotf. corr. doto HlPO z' extinc. ond flotf. corr. data 








0.72 0.74 0.76 0.78 0.80 0.82 0.84 

Time JHD - 56566 


0.72 0.74 0.76 0.78 0.80 0.82 0.84 

Time JHD - 56566 


Fig 1 Extinction and flat field corrected absolute photometry of HIPO’s B (left) and z' (right) 
channels before removal of correlated noise: the transit is clearly visible with high confidence - 
instant evidence that SOFIA is indeed a great platform for photometric observations of this kind. 
For this first iteration fit (blue) the “starspot” region (red) was excluded. Green: an alternative fit 
with fixed limb darkening. 

This noise model was computed independently from the later MCMC analysis (see section 
13.51) and no marginalization has been done over the instrument model correlation slope terms. 
We choose this approach because the lack of baseline would have caused convergence issues and a 
high risk of running into degeneracies between the noise model and the actual transit depth. For the 
analysis of the cycle 2 data with more baseline and house keeping data we are working on various 
integrated approaches using either Gaussian processes or a marginalization over the parameter 
regression. For this analysis we estimate the maximum errors introduced very conservatively by 
calculating the difference between the standard deviations before and after correction (see Table [5]) 
and using it to derive the outer error bars in Figure |71 
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Fig 2 Time series of observational parameters for the B (blue) and z' (red) ehannel. From top left 
to bottom right (for means and details about the deeorrelation see Table [B : airmass, X position 
of target star on deteetor, Y position, FWHM of PSF, altitude, maeh number, density, foeus, PSF 
eeeentrieity, oversean level, sky baekground, magnitude of eomparison stars, X position of eom- 
parison star 1, Y position of eomparison star 1, X position of eomparison star 2 and Y position of 
eomparison star 2. The position the eomparison star has a different trend than the target sinee it 
rotates around the Foresight, diseontinuities in some parameters are eause by the unavoidable line 
of sight rewinds. 
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As an alternative we performed a prineipal eomponent analysis (PCA) on this timeseries of ob¬ 
servational parameters and used the prineipal eomponents as veetors Pi{t) in our modeling of the 
residual noise. For a PCA the veetors of the observational and environmental parameters (see Fig- 
ureO that deseribe the temporal evolution of the observational environment (e.g., Maeh number or 
flight altitude) or instrument behavior (PSF shape, deteetor position) are orthogonally transformed 
to eonvert this eorrelated set of variables into a set of values that are linearly uneorrelated - the 
prineipal eomponents. Advantages of using a PCA to disentangle the observational parameters are 
that it breaks degeneraeies between parameters and reduees the number of fitting parameters. A 
disadvantage is loss of physieal insight, sinee prineipal eomponents are not always direetly eon- 
neetable to the original parameters. 

The linear eombination of observational parameters and PCA approaehes were similarly effee- 
tive in modeling the residual noise. The linear eombination approaeh produeed a slightly better 
rms of the final light eurve so we used this method with the original veetors for the rest of our 
analysis. 
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Table 1. Correlation and fit coefficients of the observational parameters with the residuals 


Parameter pi 

Correl. (B) 

fit coeff. Ci (B) 

mean (B) 

Correl. {z') 

fit coeff. Ci {z') 

mean {z') 

Airmass 

-0.30 

-0.003 

1.91 

0.01 

-0.007 

1.847 

X position HD 189 

-0.04 

0.017 

288.66 px 

0.02 

0.15 

310.24 px 

Y position HD 189 

0.30 

0.048 

833.130 px 

-0.12 

-0.155 

839.67 px 

FWHM 

0.06 

-0.009 

4.259 px 

0.03 

-0.002 

3.968 px 

Altitude 

-0.11 

-0.122 

38995.1 ft 

0.03 

0.092 

38994.7 ft 

Mach number 

0.03 

0.0155 

0.847 

-0.11 

-0.008 

0.848 

Density 

0.17 

-0.036 

0.305 kg/m^ 

-0.20 

-0.04 

0.306 kg/m? 

TA Focus 

-0.12 

-0.009 

852.61 u 

0.46 

0.006 

846.98 u 

Image eccentricity 

-0.09 

-0.005 

0.111 

-0.26 

-0.0002 

0.049 

Position angle 

-0.24 

0.0005 

-32.87 deg 

-0.06 

-0.00003 

-56.683 deg 

Overscan 

0.30 

0.714 

1011.75 ADU 

0.06 

0.009 

1419.01 ADU 




Table 1—Continued 


Parameter pi 

Correl. (B) 

fit coeff. Ci (B) 

mean (B) 

Correl. {z') 

fit coeff. Ci {z') 

mean {z') 

Sky Brightness 

-0.36 

-0.028 

5.45 ADU 

0.15 

-0.00006 

62.68 ADU 

Comparison stars mag 

0.03 

-0.00007 

12.66 

-0.00 

-0.081 

10.43 

X pos, comparison 1 

-0.23 

0.0025 

202.637 pix 

0.13 

-0.168 

216.338 pix 

Y pos, comparison 1 

0.10 

-0.018 

540.079 pix 

-0.09 

0.169 

503.724 pix 

X pos, comparison 2 

-0.29 

-0.0026 

599.833 pix 

0.12 

0.024 

607.72 pix 

Y pos, comparison 1 

-0.25 

0.016 

411.21 pix 

0.00 

0.085 

391.184 pix 





3.4 Analysis of a systematic effect during the transit baseline 


As mentioned before, we initially exeluded a number of datapoints showing a systematic offset 
(that at first was identified as a potential starspot feature, red data points in Figure [B from our 
analysis, in order not to underestimate the resulting transit depth. In order to determine the source 
of this signature we looked for any potential dependence of the feature on our observational pa¬ 
rameters. The density and overscan variations appeared to have a correlated behavior to the feature 
(see Figures [3] and H]) and also occur during an extremal value of the focus parameter (see Figure 
O. To test this, we interpolated the noise models from the fits to the masked-out data onto the data 
in the region of the potential starspot feature (light blue models in Figures 0] and |3]). As an alterna¬ 
tive approach we applied an identical fit of a linear combination of the observational parameters to 
the unmasked residuals and then subtracted this noise model from the data. Using this method we 
were able to almost fully correct for the feature (green models in Figures [3] and ID). We can show 
that both models - the one based on the full light curve as well as the model extrapolated from the 
fit to the masked data - can correct for all of the spot in blue and most of it in red (see Figures IH 
and IB- A small signature in red remains, but the morphology still correlates with the overscan, 
density and focus parameters. 

Considering the correlation with observational parameters and a wavelength dependence {zt 
larger than B) that is opposite than expected from a starspot occultation. We therefore conclude that 
this light curve feature is most likely correlated with a non-linear dependence of the photometry 
with the observational parameters and is not caused by an astrophysical signal (as for example a 
starspot). 
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Fig 3 Top: residuals of the first order HlPO B lighteurve - the eorrelated feature (red) was exeluded 
from first order noise models - a linear eombination of observational parameters (green) or their 
prineiple eomponents (light blue). Timeseries of density (eenter top) and oversean (eenter bottom) 
illustrates the eonneetion between the light eurve feature and these observational parameters. Bot¬ 
tom: Final full residuals for HlPO B - in this ease correeted by the green model extrapolated from 
the out-of-spot data. The potential spot signature (red) can almost completely be resolved by the 
noise model. 
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Fig 4 Same as Figure [3] but for HlPO z'. 
3.5 MCMC Transit Fit 


To determine the uneertainty in the light eurve parameters properly, we analyzed the eorreeted 


transit lightcurves with the IDL Transit Analysis Paekage (TAP) (l4ll) . The TAP software utilizes 


Markov Chain Monte Carlo (MCMC) techniques to fit transit light curves using the (l39t) model. 


TAP uses EXOFAST (l42t) to calculate the models. The package incorporates a wavelet-based 


likelihood function developed by (l43h which allows the MCMC algorithm to assess parameter 
uncertainties more robustly than classic methods by parameterizing uncorrelated ‘white’ and 
correlated ‘red’ noise. However, we did not use the TAP red noise modeling in our reduction, but 
instead tried to model the correlated noise using our observational parameters (see discussion in 
Section [33l) . 
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The planet-to-star radius ratio {Rp/R^), the time of mid-transit (Tq) and the above-mentioned 
white noise and airmass parameters were the only free parameters in our fit. We ran two iterations, 


one with the limb darkening parameters fixed to values from (1441) and the other with the limb 
darkening parameters as fitted parameters. 

(1451) showed that fixing the limb darkening eoeffieients to theoretieal values results in inaccu¬ 
racies in the fitted planetary radius of more than 1-10%. They also demonstrated how the presence 
of stellar spots can be responsible for this disagreement between the limb darkening tables and the 
fitted limb darkening coefficients, but showed that this does not affect the precision of the planet 
radius determination. Our results in comparison to previous observations confirm this finding. The 


recent analysis of (1461) confirms that fixing the limb darkening coefficients can lead to significant 


biases up to a few percent and also conclude that, in this case, the best approach is to let them be 
free in the fitting procedure. 

The scaled semi-major axis (a/i?*), the inclination (i)), the orbital period (P), the eccentricity 
(e) and the argument of periastron (ce) were fixed to previously published values (e.g.. Table O 
because these parameters are usually difficult to constrain using only a single transit for each target 
and are very accurately determined for HD 189733 b. For each transit light curve, we ran 5 MCMC 
chains of 10® steps each. 


3.5.1 Analysis of the Full Lightcurve 

We showed that using a noise model, fitted as a linear combination of the observational vectors, 
enabled us to correct for most of the noise signature that looked like a potential starspot signature 
(see 13.41) . Therefore we had to assume that it was most likely an instrumental effect and that we 
can also test our method on the whole data-set, including the potential starspot signature. For this 
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analysis we used the same method as deseribed above, but fitting our noise model and the transit 
lightcurve to all available data. Like the results shown in Figures [3] and |4] the full model was able 
to eorreet for the signature in the B ehannel but only partly in the z' ehannel. The fits to the full red 
lighteurve were elearly mislead by the above baseline values of this noisy signature (see e.g. Figure 
HI) leading to a mueh shallower estimate of the z' transit depth (see results in Table H]). Without 
more baseline data and eorresponding housekeeping parameters it is hard to pin down the cause for 
this signature: since it strongly correlates with observational parameters (see Table [Hand Figures 
lUandlll) we suspect non-linear effects that were left uncorrected from our linear noise model. In 
particular the TA focus value with a very high correlation value of 0.46 in z' (see Tabled]) and an 
extreme value during the time of the signature (see Figure [2l center) is known to have non-linear 
effects on the photometry. Therefore we chose to prefer the results of the lightcurves with these 
data excluded for our futher analysis. 

4 Results 

Our fit results for the corrected HIPO B lightcurves (see Figure [5|) are shown in Tabled and our fit 
results for the corrected HIPO z' lightcurves (see Figure are shown in Table (3] We find transit 
depths in B and z' that are consistent with previous observations (see Figure |7]) with HSTs ACS 
and STIS instruments and confirm the slope in the optical most likely due to some form of Rayleigh 
scattering. 

We believe the accuracy of future transit observations with SOFIA will be significantly im¬ 
proved compared with our current results (more in Section [5^ : Our observation suffered from 
a lack of out-of-transit baseline, decreasing the overall precision of the transit measurement, and 
there is also room for improvement in the understanding of the systematics caused by the air- 
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Table 2. MCMC fit results for the HIPO B lighteurve, ‘starspot’ excluded 


Parameter 

raw, fixed ED 

raw, free ED 

corr, fixed ED 

corr, free ED 

Period (d) 

2.2185741’^ 

2.218574D 

2.218574D 

2.218574D 

Inclination (°) 

85.591’^ 

85.59D 

85.591’^ 

85.591’^ 

a/R* 

8.8251’^ 

8.825D 

8.825D 

8.8251’^ 

Eccentricity 



O.'" 

0.^ 

Omega (°) 




0.^ 

j 

n 1 +U.UUU&3 

O.iJJOJ _o.00054 

0.15686^^4 

A -I CiJOT +U.0U040 

W.iJJZ/ _0,00041 

0 1 56 1 4 +0-00052 

_0.00060 

Mid Transit 

U./OOJi _o,00011 

0.78629 

A 1 +0.000082 

U. / OOJt i _o,000081 

A 7C^;48S +0.000084 

U. / 00487 -0.000084 

Einear ED 

0.8541® 

0.978 

0.8541® 

0 967 +0-024 
U.yo/ _0.049 

Quad ED 

-0.00921® 

-0 1 55 +0-047 
U.IJJ _o,030 

-0.00921® 

A -1 OQ -1-0.064 
-0.178 _o,034 

Airmass Y-int 

0 QQQ15 +0.00021 
O.VVV7J -0.00021 

0 99940 +0-00024 

-0.00023 

0 QQQ02 +0-00016 
\j.yyy\jz _o.oooi6 

0 QQQ05 +0-00019 
v.yyyvj -o.ooois 

Airmass Slope 

0 0103 +0 0036 
U.U1U7 _0.0036 

0 0097 +0-0037 
yj.yjyjyi _o.oo37 

0 0161 +0-0028 
U.UlDi _0.0027 

0 0153 +0-0029 

V.VLJJ _0.0029 

Sigma Red 

0." 

0." 

0." 

0." 

Sigma White 

0 001 145 +0-000037 

U.UUllH-J _0.000035 

0 001 137 +0-000037 
U.UUllJ/ _0.000035 

0.000879 

0.000877 


®^Value Fixed from 


m 


^Value Fixed from (l44t) 


“^Fixed in first fit in order to model red noise independently 


borne environment. Additionally, the observations were further degraded due to the low reflectivity 
(30%) of the current SOFIA tertiary mirror. Future observations with HIPO will seek to improve 
upon the first two areas, and we hope to reach a similar precision in the infrared with FLITECAM. 

4.1 Comparison to Expected Photon Noise 

We performed a calculation of the expected photon noise levels based on our observed star and 
sky signal levels. We find a signal in electrons per exposure of 5.81 x 10® in the B channel,and 
1.44 X 10^ in z'. The total sky contribution in this aperture in electrons is 1.46 x 10® for B, 
and 2.83 x 10® in z'. The total noise was computed by taking the square root sum of the shot 
noise values of the stellar signal (2410 e“ B, 3800e“ z’'), the shot noise contribution of the sky 
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Table 3. MCMC fit results for the HIPO z' lighteurve, ‘starspot’ exeluded 


Parameter 

raw, fixed LD 

raw, free LD 

eorr, fixed LD 

eorr, free LD 

Period (d) 

2.218574D 

2.218574D 

2.218574D 

2.218574D 

Inelination (°) 

85.59D 

85.59D 

85.59D 

85.59D 

a/R* 

8.825D 

8.825D 

8.825D 

8.825D 

Eoeentrieity 


0.*" 

0.*" 

0.^ 

Omega (°) 


0.*" 

0.*" 

0.^ 

Rp/R*"" 

0 1 3431 +u.uuui8 

O.lJHJl _o .00018 

A -1 C0 1 c +U.0U034 

O-rJ-JlJ -0.00046 

A -I o +u.ouoi<' 

0.19418 _0.00018 

A -| +U.00053 

O.IJJJJ _o.00049 

Mid Transit 

A 7csa.a.i +0-000033 
U./0D441 _0.000033 

A +0-000033 

0. /034dU _0.000033 

A 7cssna. +0.000031 
0. /899U4 _o. 000033 

n 7S^^9r> +0.000031 

U. /899ZU _o.000031 

Linear LD 

0.349^ 

A -1 oo +0.089 
O-r-J-J -0.082 

0.349^ 

0 219 +0-087 
O.ZIV _o,083 

Quad LD 

0.221^ 

0 41 2 +0-094 

O.H-iZ _Q ;^Q 

0.221^ 

0.285 tofo^ 

Airmass Y-int 

0 QQQQ10 +0-000057 

O.VVVVIO _o .000057 

0 QQQQ^O +0-000059 

v.yyyyjv _o.oooo 60 

0 QQQR14 +0.000050 
U.VVV07+ -0.000053 

0 QQQQ2Q +0.000056 
yj.yyyyzy _o.oooo56 

Airmass Slope 

0 0017 +0-0010 

0.00 i / _0.00099 

A 009^9 +0.00099 
V.VVZJZ _0.00099 

0 0031 1 +0.00091 
U.UU911 _0.00097 

0 00989 +0.00094 
U.UU28Z _o.oo093 

Sigma Red 

O.'i 

0.^ 

O.'i 

O.'i 

Sigma White 

0.000728 

0 000721 +0-000011 
U.UWU/Zi _o.000011 

0 000094 +0-000011 

U.UWUDV+ _o,000010 

0.000681 


^Value Fixed from (l40h 


^Value Fixed from (l44t) 


‘^0.0012 too low at z' due to eompanion eontamination 


'^Fixed in first fit in order to model red noise independently 


Table 4. MCMC fit results for the full HIPO lighteurve, systematie outlier ineluded 


Parameter'^ 

HIPOB 

HIPOz 

j 

0 15607 +0-00054 
O.iJDO/ _o.00059 

n 1 ^97^ +0.00025b 
V.iJZ/J _o.00022 

Mid Transit 

0.786505 1S:»S33 

0 78SS4S +0.00003 

U. /89949 _o.oo003 

Linear LD 

0.969 1^22 

0 033 +0-043 

V.VJJ _0.024 

Quad LD 

1 oo +0.063 
-.IJJ -0.033 

0 514 +0-029 

W.J1+ _0.05 

Airmass Y-int 

0 QQQ03 +0-00021 
v.yyyvj _o.oo02l 

OQQQQRl +0.00006 
v.yyyyol -0.00006 

Airmass Slope 

0 014 +0-003 

U.UIH- _0,003 

o 

o 

o 

4Z 

1 + 
o o 
b b 
o o 

Sigma White 

A AA1 +0.00003 
U.UUi _o.00003 

0.000768 


^All other values fixes from 


m 


'^0.0012 too low at z' due to eompanion eontamination 
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Fig 5 Final lightcurve in HIPO B, fit (blue, with free limb darkening). Triangles are averaged over 
30 datapoints. For derived preeision (also for the binned data) see Table [51 



HIPO B final data/fit 



baekground (380 e“ B, 1680 e~ z'), and a read noise (6e“ per pixel) of 975 e~ on both sides. 

This gives a total noise of 2600 e~ in B and 4300 e~ in z', resulting in a final photon noise 
limited signal-to-noise ratio (SNR) of 2200 per 7 seeond integration in B, and 3400 per 3 second 
integration in z'. This is very close to the theoretically derived values of a SNR of 2010 in B and 
3990 in z' from the HIPO exposure time calculator spreadsheet (provided by the HIPO instrument 
PI q). Most of the difference is due to the sky contribution, which was unassessed prior to the 
observations. With our observation we were able to obtain real data on the sky background at 


SOFIA’S flight altitude, which will be used to improve future exposure time estimates. 


Adding an allowance for scintillation noise (a poorly understood contribution that may be im- 
*http;//www. sofia.usra.edu/Science/ObserversHandbook 
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HlPO z' final data/fit 
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Fig 6 Final lightcurve in HIPO z\ fit (blue, with free limb darkening). Triangles are averaged over 
87 datapoints. For derived preeision (also for the binned data) see Table [51 


portant for bright stars like HD 189733 - here we used the recipe from the aforementioned HIPO 


exposure time calculator spreadsheet that is based on the findings in (l2flt) f. the derived SNR values 
are adjusted downward to about 1900 in B and 2800 in z'. We used these numbers for comparison 
and also to model photon noise limited lightcurves to test the MCMC fitting routine. Indeed, we 
find that the limited amount of out of transit data increases the error by a factor of 1.5. The final 
photometric precision (600 ppm in B and 530 ppm in z' for the fitted parameter {Rp/R^), see 
FigurejTl corresponding to 187 ppm in B and 165 ppm in z! for the actual transit depth {Rp/R^Y 
respectively) is close to two times the expected photon noise (see error bars in Figure |7]). While 
our initial noise model (l3.3l) was able to reduce the raw standard deviation araw by about 50% in 
B and 30% in z' (see Table [5]), we are still not at the photon noise lever yet. As already discussed 


24 













Fig 7 Fit results of our HIPO-SOFIA observation (blaek: after noise eorreetion, grey: maximum 
uncertainties of noise model)in comparison to previous observations from other space-based plat¬ 


forms such as HST and Spitzer (see legend, from 12a; 1341) . The error bars are given by the TAP 


MCMC routine. Our B and z' data points reproduce the slope in the optical, most likely caused by 
Rayleigh scattering, though with an overall smaller transit (^th, with might be caused by a differ¬ 
ence in unocculted starspots between the observed epochs (l26b . Given that our observation (1) had 
very little out of transit baseline, (2) taken with an only 30% reflective tertiary and that (3) there is 
still room for improvement for the correction of systematics once they are better understood with 
further observations, we can report that SOFIA should be able to deliver space based data quality 
at least in the optical with HIPO. The goal for the future is to reach similar quality in the infrared 
with FLITECAM, to cover a wavelength region that will not be available after the end of Spitzer 
and the decommissioning of NICMOS until the start of JWST. 
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Table 5. Error analysis. Resulting preeision for both ehannels: fit parameter A{Rp/R^,) - 
theoretieal value and from EXOFAST MCMC; standard deviations for the binned data (from 
Figures [5] and [^, the final data after the slope eorreetion, and the raw lighteurves; final errors for 

the transit depth 



A{Rp/R,) 


stand, dev. 



A{Rp/R,f 

Channel 

theoretieal 

from MCMC 

^binned 

^ final 

^raw 

from MCMC 

B 

0.00051 

0.00060 

0.00046 

0.00102 

0.00157 

187 ppm 

z' 

0.00036 

0.00053 

0.00033 

0.00094 

0.00118 

165 ppm 


briefly in l3.4l we suspeet non-linearities, in partieular in highly eorrelated observational parameters 
sueh as oversean, density or TA foeus (see Table [Hand Figures [3] and [^ to be responsible for this 
residual eorrelated noise. 


5 Summary and Conclusions 

The results of our observation ean be summarized as follows: We demonstrate absolute photometry 
in the airborne environment with the HIPO instrument and show that SOFIA is a quasi-spaee based 
platform for photometry in the optieal. We have a reasonable understanding of the systematie 
noise of photometrie observation in SOFIA’S airborne environment, but there is still room for 
improvement. The biggest improvement would be a larger amount of baseline data to use for 
deriving the systematie error eorreetions. 


5.1 HD 189733 b 


A study of (l34t) provided the transmission speetrum aeross the entire visible and infrared range in 


a eonsistent analysis, ineluding the general system parameters and stellar limb darkening. Their 
resulting speetrum does not eonfirm the predieted eloud-free atmospheres of highly-irradiated hot 


Jupiters. Instead they show a Rayleigh seattering slope over the whole visible and near-infrared 
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range, with only narrow Na and K lines deteeted. (l26t) . however, argued that the slope between 


the visible and near-IR shown in Figure 9 ean partly be eaused by unoeeulted starspots and not 
Rayleigh seattering alone. 

We were able to eonfirm this previous reported slope (l34h . that ean be explained by a haze of 


eondensate grains, in the upper atmosphere of HD 189733 b (see Figure|7]). The slight overall offset 
between the SOFIA and HST observations, however, eould indeed be explained by a differenee in 
starspot eoverage between the epoehs and eonfirm the (l26t) interpretation that the slope is eaused by 


unoeeulted starspot eombined with Rayleigh seattering from moleeular hydrogen only instead of 
dust. Furthermore a differenee in spot eoverage ean also explain the differenee between the fitted 


free limb darkening values from t 


depths in the eorresponding fits (145 


le fixed theoretieal values and the divergent resulting transit 


4 ^. 


5.2 Comparison to other platforms 

Our transit observation of HD 189733 b with SOFIA-FLIPO eover some of the bands that were 
observed with other platforms {HST and Spitzer). This enables us to eompare SOFIA with the 
best existing observatories, paving the way for future observations of a wide range of exoplanet 
atmospheres with FLIPO on SOFIA. 

The eomparison of the sensitivity of our observation (187 ppm in B and 165 ppm in z' for 
the transit depth {Rp/Rff) with similar observations with HST and Spitzer (see Figure |7]) show 


that SOFIA is a eompetitive platform in this fie 


d, inpartieular if viewed next to to the very first 


exoplanet observations with these platforms (Ul; l47|; 1481). Comparing these early proof of eoneept 


observations with this very first SOFIA exoplanet observation illustrates this. In detail the room 


for improvement is threefold: (1) as deseribed earlier the observatin was limited to less than a 
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quarter of an hour out of transit baseline, about a faetor 6-8 less data then aetually expeeted ( 2 ) 
only about 30% of the light was passed to HIPO by the refieetive tertiary, whieh provides room 
for about 3 times more signal (3) lessons learned from this and upeoming exoplanet observations 
(Pis: Angerhausen/02-0053, 03-0052, Dreyer/02-0084, Huber/03-0042 and Swain/03-0037) will 
further improvement our understanding for the eorreetion of systematie noise. 


5.3 Future prospects with SOFIA 

The main goal for the future has to be to reaeh similar quality in the infrared with FLITECAM. 
This will eover a wavelength region that will not be available after the end of Spitzer and the 
deeommissioning of NICMOS until the start of JWST. In partieular the planned observation to 
demonstrate FLITECAM’s speetroseopie transit mode (PI: Swain/03-0037) will hopefully pave 
the way for this. 

A fully refieetive tertiary and speeiality filters for HIPO (e.g. speeial filters to avoid O 3 and 
H 2 O, or to analyse narrow band K or Na) will further improve SOFIA’S eapabilities to eharaeterize 
exoplanets. In eontrast to spaee-teleseopes it is possible and highly desired to update SOFIA with 
state-of-the-art instrumentation. With regard to upeoming ealls for next generation instruments, we 
support the idea of equipping SOFIA with a dedieated 2”'^ generation exoplanet instrument (e.g.. 


NIMBUS l49t). Another important point is to ensure better baseline eoverage for future observation, 
e.g. by allowing speeial deployments for these kind of time eritieal observations. 

In the eontext of other platforms SOFIA has obvious synergies with ground-based teleseopes: 
SOFIA is ‘filling the gaps’ between bands that are unobservable from the ground due to mostly 
H 2 O absorption. These are erueial bandpasses to eharaeterize important moleeules (H 2 O, CH4, 
CO 2 , PAHs ete) in exoplanetary atmospheres. In partieular after the rejeetion of eharaeterization 
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proposals for dedicated space observatories (like EChO or FINESSE), SOFIA is - until the start of 
JWST - the only quasi space based platform >1.7 micron {HST WFC3 limit). 
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List of Figures 

1 Extinetion and flat field eorreeted absolute photometry of HIPO’s B (left) and z' 
(right) ehannels before removal of eorrelated noise: the transit is elearly visible 
with high eonfidenee - instant evidenee that SOFIA is indeed a great platform for 
photometrie observations of this kind. For this first iteration fit (blue) the “starspot” 
region (red) was exeluded. Green: an alternative fit with fixed limb darkening. 

2 Time series of observational parameters for the B (blue) and z' (red) ehannel. From 
top left to bottom right (for means and details about the deeorrelation see Table [B : 
airmass, X position of target star on deteetor, Y position, FWHM of PSF, altitude, 
maeh number, density, foeus, PSF eeeentrieity, oversean level, sky baekground, 
magnitude of eomparison stars, X position of eomparison star 1, Y position of 
eomparison star 1, X position of eomparison star 2 and Y position of eomparison 
star 2. The position the eomparison star has a different trend than the target sinee 
it rotates around the boresight, diseontinuities in some parameters are eause by the 
unavoidable line of sight rewinds. 


37 


3 Top: residuals of the first order HIPO B lighteurve - the eorrelated feature (red) 
was excluded from first order noise models - a linear combination of observational 
parameters (green) or their principle components (light blue). Timeseries of density 
(center top) and overscan (center bottom) illustrates the connection between the 
light curve feature and these observational parameters. Bottom: Final full residuals 
for HIPO B - in this case corrected by the green model extrapolated from the out- 
of-spot data. The potential spot signature (red) can almost completely be resolved 
by the noise model. 

4 Same as Figure [3] but for HIPO 2 ;'. 

5 Final lighteurve in HIPO B, fit (blue, with free limb darkening). Triangles are 
averaged over 30 datapoints. For derived precision (also for the binned data) see 
Table [5l 

6 Final lighteurve in HIPO z', fit (blue, with free limb darkening). Triangles are 
averaged over 87 datapoints. For derived precision (also for the binned data) see 
Tabled 
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7 Fit results of our HIPO-SOFIA observation (black: after noise correction, grey: 


maximum uncertainties of noise model)in comparison to previous observations 


from other space-based platforms such as HST and Spitzer (see legend, from l2^; 
341) . The error bars are given by the TAP MCMC routine. Our B and z' data 
points reproduce the slope in the optical, most likely caused by Rayleigh scatter¬ 
ing, though with an overall smaller transit depth, with might be caused by a dif¬ 


ference in unocculted starspots between the observed epochs (l26t) . Given that our 
observation (1) had very little out of transit baseline, (2) taken with an only 30% 
reflective tertiary and that (3) there is still room for improvement for the correction 
of systematics once they are better understood with further observations, we can 
report that SOFIA should be able to deliver space based data quality at least in the 
optical with HIPO. The goal for the future is to reach similar quality in the infrared 
with FLITECAM, to cover a wavelength region that will not be available after the 
end of Spitzer and the decommissioning of NICMOS until the start of JWST. 
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5 Error analysis. Resulting precision for both channels: fit parameter - 
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